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Abstract 

We will obtain real and complex solutions to the Painleve IV equation through supersymmetric 
quantum mechanics. The real solutions will be classified into several hierarchies, and a similar 
procedure will be followed for the complex solutions. 

1 Introduction 

The Painleve equations can be seen as the nonlinear analogues of the classical linear equations associated 
to the well known special functions T,'?. They have been identified as the most important non-linear 
ordinary differential equations 3 . Although discovered from strictly mathematical considerations, nowa- 
days they are widely used to describe several physical phenomena [1]. In particular, the Painleve IV 
equation (Piv) is relevant in fluid mechanics, non-linear optics, and quantum gravity (SI. 

On the other hand, since its birth supersymmetric quantum mechanics (SUSY QM) catalyzed the 
study of exactly solvable Hamiltonians and gave a new insight into the algebraic structure characterizing 
these systems. Historically, the essence of SUSY QM was developed first as Darboux transformation in 
mathematical physics ^6 and as factorization method in quantum mechanics [71[S]. Moreover, through 
SUSY QM one can obtain quantum systems described by second-order polynomial Heisenberg algebras 
(PHA), whose Hamiltonians have the standard Schrodinger form and their differential ladder operators 
are of third order. It has been shown that there is a connection between these systems and solutions 
gix;a,b) oi Pjv 

The Piv solutions can be grouped into several hierarchies, according to the family of special functions 
they are related with. This classification can be easily done for the class of real solutions [9 , but it can 
be as well performed for the recently found complex solutions |10| . which is our aim here. To do that, we 
have arranged this paper as follows: in Section 2 we shall present the general framework of SUSY QM 
and PHA. In the next Section we will generate the real and complex solutions to Piv] then, in Section 
4 we will study the real solution hierarchies and we shall analyze the domain of the parameter space 
(a, b) where they are to be found. In Section 5 we do the same for the complex solution. We present our 
conclusions in Section 6. 



2 General framework of SUSY QM and PHA 

In the fc-th order SUSY QM one starts from a given solvable Hamiltonian 
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and generates a chain of first-order intertwining relations jlffll3) 



HjA+ = A+H,^u H,^,AJ=AJH„ (2) 



By plugging Eqs. (jS)) into Eq. ([2]) we obtain 



j = l,...,k. (3) 



We are interested in the final Riccati solution afe(a;,efc), which turns out to be determined either by k 
solutions ai{x,ej) of the initial Riccati equation 

a-^{x,ej)+a\{x,€j)^2[Vo{x)~€j], j = l,...,k, (5) 

or by k solutions Uj oc exp(J 0:1(2;, ej)dx) of the associated Schrodinger equation 

HoUj ^ -^Uj +Vo{x)uj ^ ejUj, j = l,...,k. (6) 

Thus, there is a pair of fc-th order operators interwining the initial Hq and final Hamiltonians Hk, 
namely, 

HkB+ = B+Ho, HoB- = B-H,, B+ = A+ . . . A+, B' ^ A^ . . . A' . (7) 

(k) 

The normalized eigenfunctions V'n of Hk, associated to the eigenvalues En, and the k additional eigen- 
states V'ej ' associated to the eigenvalues ej which are annihilated by B^ {j ~ 1, . . . , fc), are given by PIT4] : 

,,.(k) _ B+i^n Wiui,...,Uj.i,u,+i,...,Uk) , . 

Note that, in this formalism the obvious restriction ej < Eq = 1/2 naturally arises since we want to avoid 
singularities in Vk{x). 

On the other hand, a m-th order PHA is a deformation of the Heisenberg-Weyl algebra of kind [TIHTB] : 

[i/,i±]=±i±, [L-,L+]=Qm+liH + l)-Q^m+liH)^PrniH), (9) 
Q,n+l{H) = L+L- = (i/ - £1) . . . (i/ - £ra+l) , (10) 

where Pm{x) is a polynomial of order m in x and £i are the zeros of Q„i+i{H), which correspond to the 
energies associated to the extremal states of H. 

Now, in the differential representation of the second-order PHA (m = 2), i+ is a third-order differ- 
ential ladder operator, chosen by simplicity as |17l : 

L+ L+L+, L+ = -I5 + /(^)] , £+ ^ 1 + g(^)£ + h{x)]. (11) 

These operators satisfy the following relationships: 

HL+ = L+{H, + 1), H,L+ = L+H [H,L+]=L+, (12) 

i/a being an auxiliary Schrodinger Hamiltonian. Using the standard first and second-order SUSY QM 
gives rise to 

/ 2 

f = x + g, h = -x^ + |- - Y - 2xg -I- a, (13) 

^-y-Y + T + ^-^+^-i (14) 

5" = ^ + |,9' + 4V + 2(.x2-«)g + ^. (15) 
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The last is the Painleve IV equation [Piv] with parameters 

a = £2+£3-2£i~l, b = -2{£2 - £3? ■ (16) 
If the £i, i = 1, 2, 3 are real, we will obtain real parameters a, h for Eq. (|15|) . 

3 Real and complex solutions to Pjv with real parameters 

It is well known that the first-order SUSY partner Hamiltonians of the harmonic oscillator are naturally 
described by second-order PHA, which are connected with Piv- Furthermore, there is a theorem stating 
the conditions for the hermitian higher-order SUSY partners Hamiltonians of the harmonic oscillator to 
have this kind of algebras (see [9]). The main requirement is that the k Schrodinger seed solutions have 
to be connected in the way 

Uj^{a~y~^ui, ej=ei-(j-l), j = l,...,fc, (17) 

where a~ is the standard annihilation operator of Hq so that ui is the only free seed. 

If ui is a real solution of Eq. © without zeros, associated to a real factorization energy ei such 
that ei < £'0 = 1/2, then all Uj are also real and, consequently, the solutions to Piv are also real. On 
the other hand, if we use the formalism as in [9] with ei > Eq, we would obtain only singular SUSY 
transformations. In order to avoid this we will instead employ complex SUSY transformations. The 
simplest way to implement them is to use a complex linear combination of the two standard linearly 
independent real solutions which, up to an unessential factor, leads to the following complex solutions 
depending on a complex constant A = A -I- (A, k S M) [18]: 

u{x; e) = e-V2 [^p, (i^, i; x^) + x A (^, |; x^)] , (18) 

where iFi is the confluent hypergeometric function. The results for the real case 119' are obtained by 
making k = and expressing A = A, with e R, as 

A = A = 2i/g^. (19) 

Note that the extremal states of H^. and their corresponding energies are given by 

ek = ei-{k-l), (20) 

£1 + 1. (22) 

Recall that all the Uj satisfy Eq. (IT7|) and ui corresponds to the general solution given in Eq. ([18]). 

Hence, through this formalism we will obtain a fc-th order SUSY partner potential Vfc(x) of the 
harmonic oscillator and a Pjv solution ei), both of which can be chosen real or complex, in the way 

Vk{x)^^-{\n[Wiu,,...,Uk)]}", (23) 
gkix-e,) = -x-{\n[iP£,{x)]y. (24) 

For fc = 1, the first-order SUSY transformation and Eq. (I24p lead to what is known as one-parameter 
solutions to Piv, due to the restrictions imposed by Eq. (I16p onto the parameters a, b of Piv which make 
them both depend on ei |20) . For this reason, this family of solutions cannot be found in any point of the 
parameter space (a, 6), but only in the subspace defined by the curve {(a(ei), &(ei)) , ei € M} consistent 
with Eqs. (|16p . Then, by increasing the order of the transformation to an arbitrary integer fc, we will 



W{ui,...,uk-i) 

VV(Ml,...,Ufe) 

V'f, oc S+e-"'/2, £2 = 
^£3 S+a+Wi, ^3 = 
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expand this subspace for obtaining k different families of one-parameter solutions. This procedure is 
analogous to iterated auto-Backlund transformations [21] . Also note that by making cyclic permutations 
of the indices of the three energies £i and the corresponding extremal states of Eqs. (|20II22|) (when they 
have no nodes), we expand the solution families to three different sets, defined by 

ai = -ei+2fc-^, 6i = -2(^ei + i^ , (25) 
02 = 2ei - fc, 62 = -2fc^ (26) 
a3 = -ei-fc-^, 63 = -2 + , (27) 

where we have added an index corresponding to the extremal state given by Eqs. (|20j|22l) . Therefore we 
obtain three different solution families of Piv through Eqs. (|18ll24p . The first family includes non-singular 
real and complex solutions, while the second and third ones can give just non-singular strictly complex 
solutions, due to singularities appearing in the real case. 



4 Real solution hierarchies 

The solutions gk{x; ei) of the Painleve IV equation can be classified according to the explicit functions 
they depend on [20]. In the real case, see Eqs. (I18I24P with the condition given in Eq. (|19p . the solutions 
are expressed in terms of the confluent hypergeometric function iFi, although for specific values of the 
parameter ei they can be reduced to the error function erf(x). Moreover, for particular parameters ei 
and Ui, they simplify further to rational solutions. 

Let us remark that we are interested in non-singular SUSY partner potentials and the corresponding 
non-singular solutions of P/y. Note that the same set of real solutions to Pjv can be obtained through 
inverse scattering techniques p] (compare the solutions of [20] with those of [9]). 

4.1 Confluent hypergeometric function hierarchy 

In general, the solutions of Piv are expressed in terms of two confiuent hypergeometric functions. For 
example, let us write down the explicit formula for gi{x]ei) in terms of the parameters ei, vi (with 
ei < 1/2 and < 1 to avoid singularities): 



, , _ 2.ir (^) [3 ifi (^, §; x') - {2e, + (^, |; x^)] 

3r(i^) ^F,{'-^,l;x^)+6.,xr{'-=l^) ,F,C-^,l;x^) 
^ 3x(2e, + l)r (1^) ,F,(1^,|;.^) 

+ 3r(i^) iFi(i^,i;a;2)+6^ixr(3^) ,F,C-^, I; x^Y ^"^^ 

The explicit analytic formulas for higher-order solutions gk{x;ei) can be obtained from expression (|24|) . 
and they have a similar form as in Eq. 



4.2 Error function hierarchy 

It is interesting to analyze the possibility of reducing the explicit form of the P/y solution to the error 
function. To do that, let us fix the factorization energy in such a way that any of the two hypergeometric 
series of Eq. ([TS]) reduces to that function. This can be achieved for ei = —(2m + l)/2, with m S N. 
By defining tp^-^ (x) = y^e^ [1 + vi erf(a;)], we can write down simple expressions for gk{x, ei) for some 



4 



Solution hierarchies for the Painleve IV equation 



D. Bermiidez, D.J. Fernandez C. 






X 



Figure 1: The P/y solutions given by Eqs. 



specific parameters k and ei: 

ffi(x;-5/2) = 
g2{x;-l/2) = 

4.3 Rational hierarchy 



4[i/i +X(p^^{x)] 
2iyix + {l + 2x'^)(p^^{x)' 

4:Vi[>^i + 6xtpi,^{x)] 

^•^1 {x)[pI^ [x) 



2vixip„^ (x) - 2v1] ' 



(29) 
(30) 



Now, let us look for the restrictions needed to reduce the explicit form of Eq. (|24|) to non-singular rational 
solutions. To achieve this, once again the factorization energy ei has to be a negative half-integer, but 
depending on the ei taken, just one of the two hypergeometric functions is reduced to a polynomial. Thus, 
we need to choose additionally the parameter i^i = or i^i — cx) to keep the appropriate hypergeometric 
function. However, ui have a zero at a; = when vi — >■ oo, which will produce one singularity for 
the corresponding Pjv solution. Hence, we should make j^i = and ei = —(4m + l)/2 with to € N. 
Departing from Schrodinger solutions with these vi, ei we get some explicit expressions for the gkix', ei) 
of the rational hierarchy: 



gi(x;-5/2) = - 
32(2;; -5/2) = - 
53(2:; -5/2) = - 



Ax 



2x^' 

Ax 



l&x^ 



3 + 4a;4 ' 
12(3x + 4a;3 +4x5) 
3 + 4a;4 ^ 9 + ISx^ - I2x'^ + Sx^ ' 



l-|-2a;2 
16a;3 



(31) 
(32) 
(33) 



which are plotted in Fig. [TJ 



4.4 First kind modified Bessel function hierarchy 



Another interesting case associated to a special function arises for ei = —to, to € N, which leads to the 
modified Bessel function of first kind. We write down an example of one solution belonging to such a 
hierarchy: 



51 (a;; 0) 













X 







(34) 
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Figure 2: Parameter space for real P/y solutions. The lines represent solutions of the confluent hyper- 
geometric function hierarchy, the black dots of the error function hierarchy, and the white dots of the 
rational and error function hierarchies. 



5 Complex solution hierarchies 

Let us study the complex solutions subspace, i.e., we use the complex linear combination of Eq. (|18p 
and the associated Pjv solution of Eq. ([24| . This allows the use of seeds ui with ei > 1/2 but without 
producing singularities. Moreover, the complex case is richer than the real one, since all three extremal 
states of Eqs. lead to non-singular complex Pjv solution families. 



5.1 Confluent hypergeometric hierarchy 

As in the real case, in general the solutions of Pjv are expressed in terms of two confluent hypergeometric 
functions. In particular, the explicit formula for the first family gi{x; ei) in terms of the parameters ei, A 
is given by 

, A [3 I; x') - (2ei + ?>)x\F^ (^, |; x^)] 

^ -3x(26i + l)iPi(l^,|;:r^) 

Once again, for all families the explicit analytic formulas for the higher-order solutions gk{x\t\) can be 
obtained through the formula (|24)) . 



5.2 Error function hierarchy 

If we choose the parameter ei — — (2m-Fl)/2 with to G N, as in the real case, we obtain the error function 
hierarchy. In terms of the auxiliary function 0a = e^ [4-1- A7r^/^erf(a;)], a solution from the third family 
is written as: 

, 4A-h4a;(/.A(a;) 

= 2A.+ (1 + 2.^)0a(.) - 

5.3 Imaginary error function hierarchy 

Different to the real case, now we can use t\ > 1/2, giving place to more solution families. This is clear 
by comparing the real and complex parameter spaces of solutions from Fig. [5] and Fig. [31 By defining 
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a 



Figure 3: Parameter space for complex solution hierarchies. The lines correspond to the confluent 
hypergeometric function, the black dots to the error function or the imaginary error function, and the 
white dots to the first kind modified Bessel function. 




Figure 4: Real (solid curve) and imaginary (dashed curve) parts of a complex solution to P/y. The plot 
corresponds to fc = 1, ei = 0, A 0, and k = 1. 



a new auxiliary function 0ja = e ^ [4 + A7r^/^erfi(2:)], where erfi(a;) is the imaginary error function^ we 
can write down an explicit solution from the third family 



51 (a:^; 5/2) 



4A(1 - x^) + 2x{-i + 2x^)(t>iA{x) 
2Ax+(l-2a;2)0jA(a:) ' 



(37) 



5.4 First kind modified Bessel function hierarchy 



Let us write down an example of the solution of this hierarchy for A = 0, k = 1,A = i, i.e., ui is a purely 
imaginary linear combination of the two standard real solutions associated to ei = 0: 



(I) [^1 {4) - i-k (4)] + 2-r (I) [/_3 (4) - /. ) 



r(l)/-i(f) + 2^r(t)/.(f) 

Its real and imaginary parts are plotted in Fig. 21 



(38) 
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6 Conclusions 

In this paper we have discussed a general method to obtain real and complex solutions of Painleve IV 
equation by using SUSY QM, which is closely related to the factorization method. Through this scheme 
we have shown that real factorization energies can be used to obtain Pjv solutions with real parameters 
a, b. We have shown the existence of more solutions in the complex case than in the real one by studying 
in detail the parameter space (a, b). 

We have classified the solutions into hierarchies arising both in the real and in the complex cases. 
Both classifications became very similar, except for a hierarchy which cannot be obtained in the real case. 
A further study of the Painleve IV equation with complex parameters is currently under study. 
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